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Abstract 

This  paper  addresses  an  innovative  method  for 
passive  sonar  signal  processing  where  it  is  required 
to  suppress  a  field  of  moving  acoustic  interferes 
while  simultaneously  enhancing  the  signal  from  a 
weak  moving  source.  Motivated  by  the  Space-Time 
Adaptive  Processing  (STAP)  technology  for  simi¬ 
lar  radar  applications,  we  propose  a  solution  to  the 
passive  broadband  sonar  problem  using  Wideband 
Space-Time  Adaptive  processing  ( WB-STAP )  in  the 
space-frequency  domain. 

The  overall  objective  is  to  develop  models  of  the 
source/receiver  dynamics  directly  into  the  problem 
formulation  and  analyze  time  varying  nature  of  the 
data  vector  and  the  covariance  matrix  by  operating 
in  the  spatio-frequency  domain.  This  allows  STAP- 
like  formulation  for  the  broadband  passive  sonar  in 
the  space-frequency  domain.  In  particular  by  em¬ 
ploying  subaperture  techniques  and  exploiting  the 
underlying  model  using  eigen-subspace  methods  it 
is  possible  to  obtain  optimum  adaptive  weight  vec¬ 
tors  that  require  significantly  less  number  of  data 
samples  compared  to  conventional  methods.  In  ad¬ 
dition,  this  formulation  allows  simultaneous  track¬ 
ing  of  all  targets  present  in  the  field  of  view  by 
jointly  estimating  their  arrival  angles  and  doppler 
parameters. 

1.  INTRODUCTION 

One  of  the  fundamental  problems  in  passive  sonar 
array  signal  processing  is  the  ability  to  suppress  a 


field  of  moving  acoustic  interferes  while  simulta¬ 
neously  enhancing  the  signal  from  a  desired  mov¬ 
ing  source  that  may  be  weaker  in  strength  com¬ 
pared  to  the  interfering  signals.  The  field  of  view  of 
the  passive  array  may  have  several  moving  acoustic 
sources  located  along  directions  {6k}  and  moving 
with  different  velocities  {vk}.  Further  the  acous¬ 
tic  field  so  generated  may  be  broadband  in  nature 
so  that  the  problem  reduces  to  essentially  “motion 
compensation”  and  coherent  integration  for  passive 
broadband  signals  in  presence  of  strong  interferers. 
Moving  sources  generate  equivalent  doppler  com¬ 
ponents  {(Jdk  }  so  that  the  unknown  parameter  set 
becomes  (6k,  c 0dk),  k  =  1,  2,  •  •  • ,  K,  where  K 
represents  the  total  (unknown)  number  of  signals 
present  in  the  field  of  view  and  6k  and  uodk  rep¬ 
resent  the  direction  of  arrival  (DOA)  and  doppler 
parameter  associated  with  the  kth  source. 

Similar  scenarios  exist  in  radar  where  a  set  of 
pulses  interrogate  a  target  that  is  moving  along  the 
direction  6k  and  extending  a  doppler  velocity  ujd 
at  the  array  [1,2].  In  that  case  Space-Time  Adap¬ 
tive  processing  (STAP)  has  been  successfully  ap¬ 
plied  to  suppress  the  dominant  clutter  and  interfer¬ 
ers  and  enhance  the  weak  signal  that  is  buried  in 
noise.  Adaptive  beamforming  techniques  require 
the  knowledge  of  the  extended  space-time  covari¬ 
ance  matrix,  and  as  a  result  sample  support  prob- 
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lem  arises  from  the  nonstationary  nature  of  clutter 
and  interferences  present  in  the  data  samples  that 
are  used  to  estimate  the  sample  covariance  matrix. 

To  address  this  problem  in  the  radar  context,  a 
variety  of  subspace  based  “sample  data  enhancing 
techniques”  are  available  [3,4].  To  be  specific,  other 
available  degrees-of-freedom,  such  as 

(i)  subspace  based  eigencanceler  methods 

(ii)  subarray- subaperture  spatio-temporal  smooth¬ 
ing  methods 

(iii)  convex  projection  approach  and  relaxed  nonex¬ 
panding  operators  to  project  the  covariance  matrix 
onto  the  intersection  of  desired  convex  sets  and 

(iv)  optimum  transmit-receiver  (OTR)  waveform  de¬ 
sign  methods  maximizing  output  signal  to  interfer¬ 
ence  and  noise  ratio  (SINR) 

can  be  incorporated  into  the  adaptive  weight  vector 
design  to  reduce  the  number  of  data  samples  re¬ 
quired  for  its  computation.  However  for  radar  ap¬ 
plications  STAP  works  mainly  under  the  narrow- 
band  assumption,  and  in  the  broadband  sonar  con¬ 
text  the  narrowband  methods  are  not  directly  ap¬ 
plicable.  By  reexamining  the  sonar  problem  in  the 
space-time  domain  it  is  possible  to  extent  the  STAP 
strategies  to  the  sonar  problem  using  Wideband  Space 
Time  Adaptive  processing  (WB-STAP). 

In  what  follows  we  briefly  describe  the  STAP 
problem,  and  motivated  by  that  approach  a  solu¬ 
tion  to  the  broadband  sonar  problem  is  outlined  in 
the  next  section  using  WB-STAP. 

2.  SPACE-TIME  ADAPTIVE  PROCESSING 

Space  time  adaptive  processing  (STAP)  has  emerged 
as  a  key  enabling  technology  for  airborne  intelli¬ 
gence,  surveillance  and  targeting  using  adaptive  MTI 
radar  [1].  With  its  ability  to  simultaneously  cancel 
sidelobe  jammers  and  competing  doppler- spread  clut¬ 
ter  returns,  STAP  has  the  potential  for  excellent 


detection  performance  over  much  of  the  available 
angle-doppler  spectrum.  In  normal  STAP  opera¬ 
tion,  an  estimate  of  the  “target-free”  total  interfer¬ 
ence  (clutter  plus  jammer  plus  thermal  noise)  co- 
variance  matrix  must  be  estimated  for  a  given  range- 
bin  under  test  in  order  to  calculate  the  optimum  set 
of  space-time  weight  vectors  [1,4].  In  practice,  this 
estimate  is  obtained  by  combining  the  returns  from 
range-bins  surrounding  the  test  cell  (with  an  ap¬ 
propriate  number  of  “guard  cells”  to  prevent  target 
range  sidelobe  spillover).  For  a  typical  surveillance 
radar,  a  minimum  of  several  hundred  independent 
and  identically  distributed  (i.i.d.)  range-bins  would 
be  required  if  full  degree-of-freedom  (DOF)  STAP 
was  to  be  employed. 

In  summary,  a  major  issue  in  space-time  adap¬ 
tive  processing  (STAP)  for  airborne  moving  target 
indicator  (AMTI)  radar  both  for  identification  and 
localization  is  the  sample  support  problem. 

In  this  context,  consider  the  radar  scenario  where 
the  returned  space-time  snapshot  signal  consists  of 
a  target  echo  and  interferences  such  as  jammer,  clut¬ 
ter  and  thermal  noise  given  by 

x2  =  a*a t  +  c i ,  (1) 

where  at  and  a.t  =  a (0*,  u)dt)  are  the  complex  tar¬ 
get  attenuation  factor  and  target  steering  vector  re¬ 
spectively  associated  with  the  spatial  and  doppler 
parameters  9t  and  ujcit  of  the  moving  target,  and  c i 
represents  the  total  interference  plus  noise  signal 
[1,4].  Here  £  (jMN  represents  the  concatenated 
space-time  data  vector  formed  from  the  array  out¬ 
put  vectors  corresponding  to  M  pulse  returns  that 
are  present  in  a  coherent  processing  interval  (CPI) 
with  interpulse  interval  equal  to  T.  Using  N  an¬ 
tenna  elements,  if  Xk(U)  represents  the  k- th  sensor 
output  at  t  =  ti ,  then 

x(tj)  =  [xi(ti),X2(ti),  ■  ■  ■  ,XN(ti)\  (2) 
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represents  the  array  spatial  output  at  time  t%  and 

Xi  =  [x(ti),  x(ti  +  T),  •  •  •  x(ti  +  (M-  1  )T)]r 

(3) 

represents  the  space-time  snapshot  data  vector  due 
to  M  pulses.  In  the  point-doppler  estimation  prob¬ 
lem  the  optimum  weight  vector  w  is  given  by  [4] 

W  =  R~la.t ,  (4) 

where 

R  =  E{  c,c*}  (5) 

is  the  total  interference  plus  noise  covariance  ma¬ 
trix.  In  actual  practice,  the  interference  covariance 
matrix  R  must  be  estimated  by  making  use  of  the 
returns  from  the  neighboring  range  bins  to  the  point 
of  interest.  This  is  most  commonly  accomplished 
with  the  formation  of  a  sample  covariance  estimate 
given  by  [4] 

1  n  1 

R=  =  -YnY*,  (6) 

n  “  J  n 

3= 1 

where 

Yn  =  [xi,x2,-",xn].  (7) 

Unfortunately,  for  the  above  procedure  to  work  it  is 
necessary  to  assume  stationarity  of  the  clutter  over 
all  range  bins  used  in  (6).  Moreover  since  the  size 
of  R  is  large  (MN  x  MN),  (6)  requires  a  very  large 
number  of  range  bins  (n  >  MN)  to  guarantee  a 
good  estimate  R  that  is  nonsingular,  and  that  much 
data  may  not  be  available  in  practice. 

It  has  been  shown  that  preprocessing  methods  - 
the  forward/backward  subaperture  methods  for  di¬ 
rect  data  reduction  and  convex  projection  methods 
for  preconditioning  sample  covariance  matrices  - 
that  help  to  significantly  reduce  required  data  sam¬ 
ples  in  estimating  the  sample  covariance  matrix  can 
be  employed  on  the  limited  data  set  for  superior 
performance  [3,4]. 


Motivated  by  the  STAP  methodology,  next  we 
address  the  adaptive  processing  for  wideband  pas¬ 
sive  sonar  signals. 

3.  SPACE-FREQUENCY  ADAPTIVE 
PROCESSING 

Consider  a  set  of  broadband  signals  s\(t),  52(f), 

•  •  •,  sx(t)  that  are  located  along  directions  9\,  62, 

•  •  •,  9k  with  respect  to  an  array,  and  moving  with 
effective  velocities  vi,  V2,  •  •  •,  vr  (along  the  wave- 
front).  A  set  of  N  passive  sensors  is  used  to  collect 
the  data  arriving  from  these  sources.  The  sensors 
can  be  arranged  in  various  configurations  (linear, 
circular,  rectangular  towed  arrays  etc.),  however, 
for  the  present  discussion  we  will  assume  a  uni¬ 
formly  spaced  linear  array  with  interelement  spac¬ 
ing  d.  Assume  that  the  source  located  along  direc¬ 
tion  6  is  moving  along  that  direction  with  an  effec¬ 
tive  velocity  v  away  from  the  array.  At  t  —  to, 
let  the  target  be  at  a  distance  r  along  that  direction 
with  respect  to  the  first  sensor.  Then  if  s (to)  rep¬ 
resents  the  target  waveform,  then  the  first  sensor 
output  x\ (t)  at  that  instant  is  given  by 

Xi(t0)  =  s(t0  -  r/c)  (8) 

and  since  the  ray  has  to  travel  an  additional  distance 
of  d  sin  9  to  arrive  at  the  second  sensor  we  get  the 
second  sensor  output  to  be 

X2(to)  =  s[to  —  (r  +  dsin0)/c],  (9) 

where  c  represents  the  uniform  underwater  acous¬ 
tic  velocity.  Notice  that  (9)  contains  information 
only  about  the  direction  of  arrival  9.  To  estimate 
the  velocity  and  the  associated  doppler  parameter, 
consider  the  array  output  at  time  t  =  to  +  T.  The 
target  moves  a  distance  v  T  along  the  ray  and  hence 
the  first  sensor  output  at  t  =  to  +  T  is  given  by 

xi(t0  +  T)  =  s[t0  +  T  -  (r  +  vT)/c\ 
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=  s[t0  -  r-  +  T(1  -  v/c)] 

—  s[to - 1-  Tcud\,  (10) 

c 

where  we  define 

A  V 

=  (11) 
c 

to  be  the  effective  doppler  velocity  associated  with 
the  moving  target.  From  (9)  the  second  sensor  out¬ 
put  at  t  =  t o  +  T  is  given  by 

^2(^0  +  T)  =  s[to  —  ^  +  -  dsinQ/c].  (12) 

Equations  (8)-(12)  can  be  used  to  determine  the  sig¬ 
nal  strength  arriving  from  K  moving  sources  lo¬ 
cated  along  directions  Q\,  62,  •  •  •,  9k,  and  moving 
with  velocities  vi,  V2,  •••,v^  away  from  the  array. 
From  (8)-(9),  the  zth  sensor  output  is  given  by 

K  d 

Xi(t )  =  V  Sk(t  -{i-  1)-  sin0)  +  Wi(t),  (13) 

ti 

where  s^it)  represents  the  reference  signal  of  the 
kth  source  at  the  first  sensor.  Similarly  using  (8)- 
(12),  we  obtain  the  output  at  time  t  +  nT  to  be 

k  d 

Xi(t  +  nT)  =  sk{t  ~  (i  ~  l)-sin0fe 
k= 1  c 

+  nTujd)  +Wi,n{t),  (14) 

i  =  1,  2, ...,  M,  n  =  1,2,  ...,7V. 

Here  n{t)  represents  the  noise  and  other  rever¬ 
beration  components  present  in  the  scene. 

Eqs  (13)-(15)  represents  the  data  available,  and 
to  make  further  progress  it  is  useful  to  represent  this 
set  of  data  in  the  frequency  domain. 

Let  Xi(uo ,  n)  represent  the  Fourier  transform  of 
Xi(t  +  nT)  in  (15).  This  gives 

K 

Xi(u,n)  =  J2Sk(L°)e~Ki^SindkejnTuJdk 

k= 1 

+  Wi,n(cj),  (15) 


where  Sk( u)  represents  the  Fourier  transform  of 
the  kth  source  output  Sk(t).  Notice  that  Eq.  (15) 
represents  the  available  space-frequency  data. 

Then  (15)  can  be  expressed  compactly  as 

K 

X(w)  =  Sk{u)A(tu,6k,udk)  +  W(u).  (16) 
k= 1 

Following  STAP  formulation  in  Sec.  2,  let 

R(cj)  =£?[X(o;)X*(a;)]  (17) 

represents  the  MN  x  MN  covariance  matrix  in  the 
space-frequency  domain.  Then  for  a  desired  target 
located  along  (0q,  Wd0)>  we  get  the  optimum  weight 
vector  at  frequency  u  to  be  [as  in  (4).] 

WH  =  R-\io)  A(LU,60,udo).  (18) 

Notice  that  direct  inversion  of  (17)  requires  at 
least  MN  samples  for  estimating  (17),  and  in  a 
nonstationary  scene  it  is  necessary  to  examine  al¬ 
ternate  procedures  to  estimate  the  adaptive  weight 
vector.  In  this  context,  if  the  signal  scene  con¬ 
sists  of  K  sources,  then  an  eigendecomposition  of 
R(o;),  as  well  as  the  adaptive  weight  vector  in  (18) 
can  be  expressed  with  r  replaced  by  K  there. 

For  example,  a  14  sensor  array  with  16  pulses 
used  in  Fig.  la  results  in  a  covariance  matrix  of  di¬ 
mension  14  x  16  =  224,  and  traditionally  atleast 
double  the  number  of  samples  are  required  for  3dB 
performance  [4].  In  Fig.  lb,  the  signal  is  buried 
in  40  dB  clutter  to  noise  ratio  (CNR)  with  —  3  dB 
SNR  and  it  is  detected  with  an  SINR  improvement 
factor  of  8  dB  using  five  samples  —  compared  to 
2  x  224  =  448  samples.  Fig.  2  shows  the  extension 
of  these  ideas  to  joint  Doppler- frequency  and  an¬ 
gle  of  arrival  estimation  for  the  broadband  passive 
sonar  problem  in  the  space-frequency  domain  for 
a  five  source  scene  with  10  dB  SNR  and  using  six 
sensors.  The  data  structure  in  (16)-(17)  is  exploited 
here  by  utilizing  the  eigensub  space  methods  in  the 
frequency  domain. 
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4.  CONCLUSIONS 


(a)  htp:  N=  14,  M=  16,  #  of  Samples  =  5,  SNR=  -3  dB,  CNR=  40  dB 


The  overall  objective  is  to  develop  models  of  the 
source/receiver  dynamics  directly  into  the  problem 
formulation  and  analyze  time  varying  nature  of  the 
data  vector  and  the  covariance  matrix  by  operat¬ 
ing  in  the  spatio-frequency  domain.  This  allows 
STAP-like  formulation  for  the  broadband  passive 
sonar  in  the  space-frequency  domain  (WB-STAP). 
In  particular  by  employing  subaperture  techniques 
and  exploiting  the  underlying  model  using  eigen- 
subspace  methods  it  is  possible  to  obtain  optimum 
adaptive  weight  vectors  that  require  significantly 
less  number  of  data  samples  compared  to  conven¬ 
tional  methods.  In  addition,  this  formulation  allows 
simultaneous  tracking  of  all  targets  present  in  the 
field  of  view  by  jointly  estimating  their  arrival  an¬ 
gles  and  doppler  parameters. 
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Range  Bin  Index  200  -1  Doppler 


(a)  Conventional  STAP 

(b)  htpsaspfb:  Ns=12,  Ms=14,  #  of  effective  samples:  90 


(b)  Sparse  Sample  STAP 


Fig.  1.  Adaptive  matched  filter  output  responses  of 
(a)  conventional  STAP  and  (b)  sparse  sample  STAP 
methods  using  five  snapshot  data  samples.  A  four¬ 
teen  element  array  with  sixteen  pulses  receives  data 
with  CNR  =  40dB.  Injected  target  with  SNR  =  -3 
dB  is  located  at  range  bin  220. 

Time-Delay  and  DOA  Estimation 


“  10CK 


Fig.  2.  Joint  doppler-frequency  and  angle  of  arrival 
estimation  of  arrivals  using  WB-STAP.  The  sources 
are  located  at  42,  151,  108,  73  and  11  degrees  with 
normalized  doppler  velocities  of  2.8,  3.9,  4.1,  6.5 
and  8.9. 
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